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We review the existing theoretical approaches to the study of the observables of hadronic 
atoms. In the relativistic perturbative approach based on the Bethc-Salpeter equation and 
Chiral Perturbation Theory we derive the general expression for the n + ir~ atom lifetime. 
The lowest-order correction to the relativistic Deser-type formula for the atom width is 
predicted to be (6.1 ± 3.1)%. 
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'. The study of the pion-pion scattering which forms one of the basic building blocks 

in the hierarchy of strong interaction processes, enables one to gain a deeper insight in 
the nature of strong interactions. According to common belief, low-energy interactions of 
pions are described within the Chiral Perturbation Theory (ChPT) II which exploits a full 
content of global QCD symmetries. In the standard scheme with a "large" condensate, 



the S'-wave nn scattering lengths are predicted to be a[] = 0.217 and a[j — ag = 0.258 in 
units of the inverse charged pion mass @. The calculations within the Generalized ChPT 
with a small quark condensate which contain more parameters, lead to a most likely 
: value a° = 0.271. Despite a significant difference between these numbers, both results for 
the scattering length a® are compatible with the experimental value = 0.26 ± 0.05 @. 
Consequently, a precise measurement of tttt scattering lengths will be an excellent test of 
the ChPT. 

The existing experimental information on the S'-wave tttt scattering lengths is ex- 
tracted either from the data on the process tcN — > tttcN or from the study of A e4 decay 
(for the review of the recent status of tttt experiments see, e.g. @). The forthcoming 
DIRAC experiment PS212 at CERN is aimed at the high-precision measurement of the 
lifetime of the ir + 7T~ atom i. This will allow for the determination of the quantity a[] — 
with an accuracy of 5% and thus will provide a decisive probe of the predictions of ChPT. 
The experiments on the measurement of the hadronic atom observables are carried out 
also at PSI, KEK, Frascati, Uppsala. 

The experimental study of 7r + 7r~ atom characteristics provides a direct information 
about the strong tttt amplitudes near threshold since the average momenta (~ am n ) of 
constituents in the atom is much less than either of two typical scales characterizing 
strong interactions in chiral theories (~ m n and 4nF n ~ 1 GeV). The huge difference 
between the bound-state and strong scales leads to a clear-cut factorization of strong and 
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electromagnetic interactions in the observables of the pionium which is already explicit in 
the lowest-order Deser formula for the ground-state atom lifetime H 



where \l/i(0) denotes the value of the nonrelativistic Coulomb ground-state wave function 
at the origin. Note that the Coulomb interaction which is mainly responsible for the 
atom formation, and strong interactions which lead to its decay, enter separately into 
the expression of the lifetime through the Coulomb wave function and the scattering 
lengths, respectively However, in the view of the forthcoming DIRAC high-precision 
experiment it is important to evaluate corrections to this formula in order to provide 
a careful determination of the strong scattering lengths. In particular, one can readily 
observe from Eq. (|T|) that the change of the "reference mass" from the charged to the 
neutral pion mass in the definition of strong scattering lengths leads to a ~ 12% variation 
in the lifetime. In order to avoid this ambiguity one has to perform a choice of the 
"reference mass" in the ideal, isotopically symmetric world and evaluate all mass shift 
corrections which emerge due to the m^i — m^o mass difference in the pionium observables. 
All other corrections of an electromagnetic or a strong origin should also be calculated. 

Historically, the first attempts to calculate the corrections to the pionium lifetime 
within the field-theoretical approaches were focused on the evaluation of the quantities 
entering into Eq. (p]): the scattering lengths and the wave function at the origin. The 
value of a(] — o,q was calculated in various low-energy phenomenological models 
The 7T7T strong interaction potential was constructed and the corrections to the pionium 
energy and the value of the wave function at the origin due to strong interactions, vacuum 
polarization, finite size of pion were evaluated, taking into account only the discrete 
Coulomb spectrum @; The vacuum polarization effect in the pionium observables was 
studied in the quasipotential approach H; Retardation correction to the pionium lifetime 
was calculated on the basis of Bethe-Salpeter (BS) approach 0; Electromagnetic radiative 
corrections to the tctt scattering lengths and the pionium lifetime were evaluated HI. Note 
however that the consistent approach to the calculation of bound-state characteristics can 
not be confined solely to the quantities entering into the Deser formula. Instead, to this 
end a systematic field-theoretical framework is needed which is free of any double counting 
problems. In the several recent papers an attempt was made to carry out a comprehensive 
field theoretical analysis of the pionium problem. Namely, in E0H9 on the basis of the 3- 
dimensional constraint theory bound-state equations the mass shift, radiative and second- 
order strong corrections have been evaluated. In our previous papers HH01 a systematic 
perturbative approach based on the BS equation is constructed, and a complete set of the 
lowest-order corrections to the pionium lifetime is evaluated. The vacuum polarization 
effect on the pionium observables is evaluated in within the nonrelativistic QED. 

In a nonrelativistic potential approach to the problem El the expression for the life- 
time, including the corrections, is written in terms of physical scattering lengths in charged 
and neutral channels in the presence of Coulomb interactions. This is a counterpart of the 
lowest-order Deser-type factorization: the bound-state observables are expressed in terms 
of the quantities characterizing the scattering process, and these quantities differ from the 
purely strong ones by a small and calculable mass shift and electromagnetic corrections. 
In order to evaluate these corrections, one needs the tttt potential to be given explicitly. 
Assuming local and energy-independent strong potentials in the isospin symmetry limit, 
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one reproduces the nir scattering phase shifts calculated from ChPT. At the next step 
the Coulomb potential is added and the mass shift effect is included via assigning the 
physical masses to pions, the strong potential remaining unchanged. The sign of the 
largest mass sh ift correction in the scattering lengths and the pionium lifetime calculated 
in this manner li^l turns out to be opposite to that emerging from the field-theoretical cal- 
culations E30>0, provided the charged pion mass is chosen to be the "reference mass". 
Owing to the derivative character of pion couplings in the chiral Lagrangian, a possible 
reason for this discrepancy might be an explicit energy independence of strong potentials 
used in these calculations. 

To summarize, one needs an accurate theoretical estimate of the lifetime of the 
pionium which will be measured in a high precision DIRAC experiment. The underlying 
dynamics of pions should be described by the Lagrangian of ChPT in order to provide a 
reliable and unambiguous test of the predictions of ChPT. 

In our previous papers &0 we have formulated a consistent field-theoretical ap- 
proach to the calculation of the pionium observables. Below the basic ideas and assump- 
tions of our approach are given. 

The pions within our approach are described by the elementary fields in the La- 
grangian. This allows one to write down the bound-state BS equation for the metastable 
7r + 7r~ atom. The necessary link to the strong amplitudes from ChPT where the pions are 
described by the pseudoscalar quark densities, proceeds via the Deser-type factorization 
in analogy with the potential theory. Namely, we express all corrections to the Deser 
formula in terms of the on-mass-shell pion amplitudes which are further identified with 
the amplitudes from ChPT. 

In our calculations we always adopt the so-called "local" approximation which 
consists is suppressing the relative momentum dependence of the strong amplitudes. The 
origin of this approximation can be traced back to the huge difference between the bound- 
state and strong momentum scales. Since the tree-level strong amplitude for the process 
7r + 7r~ — ► 7T 7r° does not depend on the relative momentum of the pions, the corrections 
arising due to this approximation, by a simple power counting, should be suppressed by a 
product of two small factors: the fine structure constant a which is the ratio of the bound- 
state momentum and the pion mass, and the factor m vr /(47ri ? V) ~ 10 _1 being the ratio 
of two strong scales in chiral theories. For this reason we at the present stage completely 
neglect the corrections coming from using the "local" approximation. 

The pionium in our approach is described by the exact bound-state BS equation 

< X \G \P) =< X \V(P) (2) 

where x denotes the BS wave function of pionium, Go is the free Green's function of 
the 7r + 7r~ pair and V is the full BS kernel which apart from all two-particle irreducible 
diagrams includes the self-energy corrections in two outgoing charged pion legs EIH The 
square of the cm. momentum P takes the value P 2 = M 2 = M 2 — iMT where M denotes 
the "mass" of an atom, and V stands for the decay width. 

In order to perform the perturbative expansion of the bound-state observables we 
single out from the full kernel V the instantaneous Coulomb part Vc which is responsible 
for the formation of the bound state composed of tc + and 7r~. The "remainder" of the 
potential which is treated perturbatively is denoted by V — V — Vc- 

The "unperturbed" part of the full kernel Vc is chosen in the form guaranteeing the 
exact solvability of the corresponding BS equation. The solution of this equation yields the 
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relativistic Coulomb wave function ipc with the corresponding eigenvalue P* 2 = M* 2 = 
m 2 (A — a 2 ) where m n is the charged pion mass. The solution of the full BS equation is 
expressed in terms of the unperturbed solution as follows 

< x| = constx < ^|[1 + (G^(P) - Go\P*) - V'{P))G R Q}- 1 (3) 

where Gr denotes the pole-subtracted part of the exact relativistic Coulomb Green's 
function and Q is the projection operator on the subspace orthogonal to the unperturbed 
ground-state BS wave function tSO. Substituting this solution into the complete BS 
equation (Q) we arrive at the final relation 

< Vto|[l + (GoHP) - G \n - V\P))G R Q]-\G \P) - G \n - V(P))IVto >= 

which serves as the basic equation for performing the perturbative expansion of the bound- 
state observables. The explicit form of the operator GrQ is given and the perturbation 
kernel V is known up to any given order of loop expansion. The only unknown quan- 
tities entering (parametrically) into this expression are the real and imaginary parts of 
P 2 . Consequently, expanding this expression in powers of V, we arrive at the recursive 
relations which define the mass and width of the bound state in any given order. 

In order to provide the necessary Deser-type factorization in the obtained results, 
one needs a further classification of diagrams entering into the definition of V . We 
consider the following subsets of diagrams: 

1. A purely strong part, which is isotopically invariant. This part survives when electro- 
magnetic interactions are "turned off' the theory. 

2. The part which is responsible for the m n ± — m n o electromagnetic mass difference 

3. Remaining electromagnetic effects, including the exchanges of virtual photons. 

The part 1 includes strong interactions which govern the decay of a pionium. The 
part 2 makes this decay kinematically allowed. Consequently, it is natural to consider 
them together, denoting the corresponding potential as Vi 2 = V\ + V 2 - The T-matrix 
corresponding to summation of the potential V\ 2 in all orders is given by Tx 2 (P) = Vu(P) + 
Vi2(P)G (P)Ti 2 (P) . The rest of the potential is referred to as V3 = V — V\i and is 
treated perturbatively. Further, the quantity GrQ in our basic relation can be written 
as GrQ = G (P*) + 5G where SG corresponds to the exchange of the ladder of Coulomb 
photons and is also considered as a perturbation. Then in the first order in V3 and SG 
the basic relation takes the form 

= -2iM*5M- < ip c \Ti 2 \ipc > + 

+ 5M < ^ o \(G^)'G T 12 \^ c > +(5M) 2 < ^ c \(l + T 12 G Q ){G^)"\ip c > + 

+ < i, c \{5M{G 1 )' - T l2 )5GT l2 \i, c > - < ^ c \(l + T 12 G )V 3 (1 + G T 12 )\^ C > (4) 

where Go — Go(M*) and the prime stands for the differentiation with respect to M*. 

In the "local" approximation T\ 2 does not depend on the relative momenta and 
thus can be taken out of the matrix elements. Neglecting first the corrections (only the 
first line of Eq. (^) contributes), we arrive at the relativistic counterpart of Eq. (ffl) 

12MX / 2 \2M*m„ ) v ; 
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Re(*T 12 ) = 167^+^-^+^(4774; 0,0) (6) 



Im(iT 12 ) = -I67r(^-^J [l-^f) (7) 

and T(s;p,q) denote the (dimensionless) S'-wave nn strong scattering amplitudes which 
include the effect of m n ± — m^o mass difference, go is the relative momentum of the 
pair at the threshold s = 4m 2 , with the magnitude given by the relation m 2 = m 2 + % 2 - 
In order to calculate the corrections to the relativistic Deser formula, one has to 
evaluate the integrals entering into the Eq. (Q). 

1. The correction due to the shift of the bound-state pole by strong interactions is 
determined by two terms in the second line of Eq. ([|). Calculating these integrals explicitly 
and using Eqs. @ we obtain 

q Aff(i) 

6 s = - o —ET- = ~ 5 - 47 x 10_3 m -( 2a o + a l) (8) 

O -C/l 

where E\ denotes the unperturbed ground-state binding energy. 

2. The correction due to the relativistic modification of the Coulomb wave function is 
given by the value of the wave function ipc the origin which for the particular choice 
of the kernel used in the calculations is related to its nonrelativistic counterpart by 

^ c (0) = tfi(0) (1-O.381a+--0 (9) 

and the correction in the decay width which is proportional to |-?/>c(0)| 2 is twice as large. 

3. The correction due to the exchange of Coulomb photon ladders is given by the first term 
in the third line of Eq. ([|) containing 5G. Using the known expression for the relativistic 
Coulomb Green's function H'0 and evaluating this integral explicitly, we obtain 

S c = (1/2 + 2.694 - lna) AE (1) /E 1 = 3.95 x 10~ 2 m T (2a° + a 2 ) (10) 

4. In order to obtain the mass shift and radiative corrections, in the perturbation kernel 
V3 one has to take into account the residual photon exchange diagram and the self-energy 
corrections in the outgoing charged pion legs. The net result of this effect consists in the 
cancellation of the correction due to the relativistic modification of the Coulomb wave 
function and in replacing the quantity in Eq. (^) by the Coulomb pole removed 
full amplitude for the process 7r + n~ — > Using the expression of this amplitude i^, 
and isolating the purely strong isotopically symmetric amplitude with a common mass 
equal to the charged pion mass, we arrive at the following expressions for the mass shift 
and radiative corrections 



_|_(_30 + 3 /Cf -/C 2 ±0 ) (12) 



where /Cf° and li denote the low-energy constants of ChPT ErE3. 

5. The correction due to the vacuum polarization is of order 0(a 2 m n /m e ) (the contribu- 
tion due to the hadronic vacuum polarization is negligibly small). Taking into account the 
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corresponding diagram in the kernel V3 and evaluating numerically the resulting integral, 
we arrive at the following result 

3 777- 

Svac = tt: a 2 — x 0.6865 (13) 
16 m e 

which completely agrees with the zero-Coulomb piece of the result given in Ref. EDI 
The additional contributions coming from one-Coulomb and many-Coulomb pieces which 
are also considered in this paper, within our approach arise in the second-order of the 
perturbation theory. 

6. The correction to the pionium decay width due to the finite size effect is caused 
by the modification of the instantaneous Coulomb interaction by the pion loop in the 
two pion-photon vertex (the corresponding diagram is contained in the kernel V3). For 
simplicity, we approximate this vertex by a monopole parameterization. After calculating 
the corresponding matrix element explicitly, we obtain 

6 F = ^ml<r 2 > v In (i m 2 < r 2 >*) (14) 

where < r 2 >y denotes the square charge radius of the pion. 

The relativistic Deser formula with all above corrections takes the form 

1 167T /2Am w \ 1 /2/ 2 . 2 j2 . _ 

7i 1 = — [^ JL j y 1 - ( a - 4) 00 (l+S S + 6 C + 5 M + 5em + <W + 6 F ) 

where a® and a 2 , denote the nir scattering lengths in the isospin-symmetric case, with the 
charged pion mass taken to be the common mass of the pion isotriplet. 

In our calculations for the constants U we take the numerical values from Ref. 
h = -2.3 ± 2.7, l 2 = 6.0 ± 1.3, l 3 = 2.9 ± 2.4, F 4 = 4.3 ± 0.9. Also, we use the values 
^f-lCt = 1.8 ± 0.9, ^-ICf = 0.5 ± 2.2 from Ref. M. Other input parameters in 
our calculations are the S-wave tttt scattering lengths: a[j = 0.217m" 1 , = —0.041m" 1 
calculated in ChPT and the e.m. charge radius of pion < r 2 >y= 0.439 Fm 0. The 
results of our calculations are shown in Table 1. 

As we observe from Table 1, the largest correction in the decay width is caused by 
the mass shift effect. Our result for the sum of mass shift and electromagnetic radiative 
corrections generally agrees in sign and magnitude with the recent field-theoretical cal- 
culations from Ref. E-3. The sign of the mass shift effect obtained in the nonrelativistic 
scattering theory approach 111 turns out to be opposite as compared to our result, and is 
of the same order of magnitude. 

Our result for the correction due to the exchange of Coulomb photons agrees with 
the result obtained in the potential scattering theory c2l and disagrees with the result from 
Refs. MM. Numerically the largest part in this effect comes from the nonanalytic lna 
piece which is exactly the same in our approach and in the potential theory, and is absent 
in Refs. HOI ' 

Thus, we have evaluated a complete set of the lowest-order corrections to the 
pionium decay width. For a full understanding of the problem, however, the reason for 
the difference in sign in the mass shift effect in the field-theoretical and potential theories 
should be investigated in detail. Also, one should have a reliable quantitative estimate 
of the accuracy of "local" approximation including a proper treatment of arising new UV 
divergences. 
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Table 1. Corrections to the ir + ir atom decay width. 



Effect 


Value 


Correction (in %) 


Strong 


h 


-0.22 


Coulomb photon exchange 


h 


+1.55 


Mass shift 


hi 


+2.99 ±0.77 


Electromagnetic radiative 


"em 


+1.73 + 2.31 


Vacuum polarization 


hac 


+0.19 


Finite size 


h 


-0.16 


Total 


O~tot 


+6.1 + 3.1 
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